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1. INTRODUCTION, BACKGROUND 


Our interdisciplinary background and fundamentally-oriented studies of the laws governing 
multi-component chemical vapor deposition (VD), panicle deposition (PD), and their interactions, 
(ca. 72 refereed experimental and theoretical publications since 1979, supported first by NASA, 
then also by DOE-METC, DOE-PETC and AFOSR) put the Yale University HTCRE Laboratory in 
a unique position to significantly advance the 'state-of-the-art' of chemical vapor deposition (CVD) 
R&D. With NASA-Lewis RC financial support, we initiated a program in March of 1988 that has 
led to the advances described in this report (Section 2) in predicting chemical vapor transport in 
high temperature systems relevant to the fabrication of refractory ceramic coatings for turbine 
engine components. This Final Report covers our principal results and activities for the total NASA 
grant of $190,000. over the 4.67 year period: 1 March 1988-1 November 1992. Sinceour 
methods and the technical details are contained in the publications listed in Section 5.2 (9 Abstracts 
are given as Appendices (Section 6)) our emphasis here is on broad conclusions/implications 
(Section 2) and administrative data (Section 3), including personnel, talks, interactions with 
industry, and some known applications of our work. 


2. RESEARCH ACCOMPLISHMENTS 

Our principal results may be conveniently subdivided into the following generic categories: 

2.1 Role of High Activation Energy Homogeneous Reactions in Limiting CVD- 
Rates and Deposit Quality for Heated Surfaces (Rosner, Collins and Castillo, 1991, 
1992,1993) 

In the presence of interfacial kinetic barriers the onset of vapor phase chemistry in the 
immediate vicinity of hot CVD surfaces is often associated with a dramatic reduction in CVD-rates 
and film quality (cf. the early investigations of Ghoshtagore et.al, 1970). To understand when 
vapor phase reactions will be "ignited”, we carried out a high activation energy asymptotic analysis 
(Castillo and Rosner, Castillo ,1992,1993), and have begun to apply it to predict threshold 
conditions for the success of CFBL predictions (eg. transport-shifted CVD phase diagrams 
(Rosner and Collins, 1989,1990,1991). In AFOSR-supported laboratory experiments aimed at 
verifying and exploiting this theory, graduate student (NASA-Training Grant fellowship holder) 
Joshua Collins observed "vapor phase ignition" (VPI) in the deposition of TiC> 2 (s) films from 
TTIP These films were partially characterized , based on visits to NASA-Lewis. We developed 
and applied the asymptotic approach of Rosner et.al. 1990, and Collins et.al. 1991, to predict 
rational onset conditions for vapor phase reactions in the CVD of TiC> 2 (s) and closely related 
oxides, and partially validated this approach in a prototypical case (eg, TiC> 2 (s) via TTIP) using a 
CVD stagnation flow reactor of our own design (Rosner, Collins and Castillo, 1993). This sets the 
stage for future use of the theory to design optimal CVD reactors, and infer homogeneous kinetic 
parameters based on well-defined VPI-threshold CVD-rate measurements in "new" (eg, Hf-, 
and/or, Zr-containing) systems. This combined theoretical and experimental CVD work formed the 
basis of the Yale U. PhD dissertation of Joshua Collins, who was supported for 4 years by NASA 
Training Grant: 50377 . 
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2.2 Prediction of Transport-'Shifted' CVD Phase Diagrams; Theory of Deposit 

Composition 

Valuable guidance for complex CVD processes has in the past been provided by 'phase 
diagrams’ which, on purely thermophysical equilibrium grounds, display the 'fields' of stability of 
various eligible surface compounds (see, eg., Spear, 1984) in terms of reactor feed conditions {eg. 
elements ratios) and surface temperature. However, as we discovered and documented for the 
CVD of Na 2 SC> 4 (l) (turbine "hot corrosion" applications) and W(s) (halogen-cycle lamps) (Rosner 
and Nagarajan, 1985; Rosner, Nagarajan, Gokoglu and Kori, 1987; Rosner, 1988, Rosner and 
Tandon,1991), significant transport-induced 'shifts' in the boundaries of these 'fields' occur due 
to chemical element segregation ('enrichment') phenomena within the nonisothermal vapor 
boundary layers adjacent to the deposit surface.Using our asymptotic CFBL multicomponent 
transport methods (Rosner, Chen, Fryburg and Kohl, 1979) including updated algorithms for 
more accurately computing the Soret effect , in this program we compared transport-shifted and 
'purely LTCE' CVD-diagrams for a prototypical ceramic coating system of special interest to 
NASA: viz, the titanium/boron system, for the case of a coaxial filament flow reactor 
configuration. Our results show how such diagrams can be used to select CVD conditions which 
will deposit the required stoichiometry at the highest possible rate. Indeed, based on recent contacts 
with US Air Force contractors (eg. Amercom, CA) this work ican be used to guide the optimum 
choice of operating conditions of CVD reactors for synthesizing near-stoichiometric SiC(s) fibers, 
as well as TiB 2 (s) coatings. Our efficient CFBL computational methods for multicomponent (>40 
vapor species) CVD systems which employ an excess of inert "carrier" gas can now be exercized 
in many new CVD situations of practical interest. 

In the deposition of doped glass for ultimate use as an optical waveguide, combustion 
products heavily laden with Si 02 -rich submicron droplets are passed over a slowly rotating surface 
which thermophoretically captures the microdroplets. The dopant (eg., Ge 02 , P 2 O 5 , B 2 O 3 -) 
content of these microdroplets, in turn, dictates the local refractive index of the ultimate waveguide 
product, and hence, must be carefully controlled and understood theoretically. We recently 
developed a fully coupled, pseudo-single phase laminar boundary layer theory to predict deposit 
compositions as a function of the mainstream loading of silicon and the dopant (here, Ge), as well 
as substrate temperature (Park and Rosner, 1989). We have shown that it is not permissible to 
neglect the effects of silica particle mass loading or even the heat of GeC >2 condensation. As a 
byproduct of our mathematical model and calculations, we have obtained the structure of the 
dopant vapor boundary layer, and, hence, the decisive ratio of dopant vapor number density 
established at the wall to that in the equilibrium mainstream. The capability of predicting both the 
rate of glass microparticle deposition and deposit composition (eg., dopant/Si element ratio) for 
each feed condition and wall/combustion gas temperature ratio should provide valuable insight into 
the optimization and improvement of this process, and generically related PD/VD processes. 
Moreover, while illustrated in this paper for the case of a single dopant (eg., Ge), our methods are 
readily extended to multi-element ceramic deposition. In principle, such information could also be 
superimposed on our transport-shifted CVD-phase diagrams (Rosner and Collins, 1989,1990; 
discussed above). 

Using similar methods, molten salt compositions (which determine corrosion potential and 
melting point) have recently been predicted in chlorine-containing environments of interest in 
proposed stationary power plants using coal-derived fuel (Rosner and Tandon,1991). 

2.3 Role of Vapor Phase Chemistry in Determining CVD-Rates; 

Vapor phase chemical reactions can play a significant role in dictating deposition rates, and 
important properties of the resulting deposit. In a previous program, using both (stagnant) "film" 
theory and laminar boundary layer (convective-diffusion) theory, we developed and illustrated a 
tractable 'asymptotic approach’ to predict multicomponent CVD-rates and dew points in high 
temperature vapor/condensate systems (see, eg, Rosner, Nagarajan, Gokoglu and Kori, 1987). 
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We calculated and systematically compared these quantities in two 'extremes' of vapor-phase 
chemical kinetics — one in which the gas-layer is taken to be in local thermochemical equilibrium 
(LTCE) everywhere along the diffusion path; the other in which the gas-layer is assumed to be 
chemically frozen (CF) everywhere. Based on this 'asymptotic' approach, we recently illustrated 
for laminar boundary layer alkali sulfate transport (Rosner and Tandon,1991, and Gokoglu of 
NASA-Lewis RC,1988) when CVD-rates will become sensitive to vapor phase chemical kinetics, 
and when, in the absence of interface kinetic barriers, a simple "chemically frozen" vapor phase 
analysis will be adequate for engineering purposes. This line of research appears to have had an 
impact on other US Govt, contractors concerned with the hot corrosion of combustion turbine 
parts (see, eg. Luthra,1989). We anticipate that our ability to calculate CVT-rates in the LTCE- 
layer limit (also useful in halogen lamp applications, Rosner, 1988) will be useful in evaluating 
metal coatings for high pressure chemical rocket engine applications {eg., Ir(s), and Re(s)) of 
interest to NASA. 

2.4 Scavenging of Vapors by Suspended Particulate Matter Within Boundary 
Layers (Castillo and Rosner, 1989; Rosner and Tandon,1991) 

In our earlier (DOE-supported) work we experimentally examined the mass transfer 
consequences of nonequilibrium scavenging of condensible vapor (Na 2 SC> 4 ) by suspended solid 
particles (MgO) within combustion gas boundary layers (BLs) (Rosner and Liang, 1988) using an 
extension of our 'flash-evaporation' technique (Rosner, 1989; Rosner and Liang, 1986) and laser 
light scattering methods. Both nonequilibrium laminar BL-theory (Castillo and Rosner, 1989) and 
our flash-evaporation light scattering experiments revealed that the alkali mass transfer rate 
reduction factor could be correlated with the scavenging particle cloud external surface area per unit 
volume. This body of experiments, and boundary layer theory (Castillo and Rosner, 1988), 
combined with our recent work on the coupling between particle inertia and thermophoresis 
(Konstandopoulos and Rosner,1991) and current theoretical studies of vapor "scavenging" by log- 
normal populations of suspended particles (see, eg. Rosner and Tandon,1991) can now be used to 
analyze and improve recent work on aerosol-assisted CVD and particle precipitation-aided CVD 
(PPCVD) (see, eg. Shimiogaki and Komiyama, 1987 and Hurt and Allendorf,1991) for the rapid 
growth of ceramic films. 

Incidentally, in the vapor/condensate equilibrium limit, our earlier theoretical model 
(Castillo and Rosner, 1989 ) has recently been successfully used to explain the rates of frost 
formation on cryosurfaces in moist air ( Epstein, et. al. ,1991)! 


2.5 Morphological Stability of Growing Deposits; Transport and Kinetic 

Phenomena Governing CVD-Coating Surface Topography and Deposit 

Microstructure 

To gain insight into the conditions governing surface roughness, especially when 
thermophoretic (or Soret-) transport is dominant, we carried out a linear stability theory in the 
(dilute) limit of low growth Peclet number and negligible interfacial energy (Castillo, Garcia- 
Ybarra and Rosner, 1990, 1991)). In this way, we identified a new type of morphological 
instability for the commonly encountered case when the deposit thermal conductivity exceeds that 
of the growth fluid medium. 

In general, it will be necessary to consider molecular transport and chemical reactions on a 
scale much smaller than that of the overall vapor diffusion boundary layer to understand the surface 
morphology and microstructure of coatings of even prescribed stoichiometry. 

By exploiting plausible hypotheses, and making use of extensive available experimental 
data for some CVD materials (eg. pyrolytic graphite and boron nitride) we have also begun to 
develop, test, and hope to ultimately recommend, a rational general "theory" for the density (and 
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other density-related properties ) of CVD materials as a function of deposition conditions. One 
promising hypothesis is that a suitably normalized density should be a universal function of a 
surface-diffusion based Damkohler number governing the adatom's search for a "crystalline" site. 
The present program provided the opportunity to demonstrate a preliminary version of this type of 
correlation, as well as extend our capability to establish mechanism/microstructure interrelations 
(Tassopoulos and Rosner,1991). 

The abovementioned concepts are now being confronted with TiC> 2 (s) deposit roughness 
and microstructural data recently obtained by J. Collins, using the Yale-HTCRE stagnation flow 
CVD reactor. These observations will be contained in his dissertation and reported in part at the 
1993 AIChE annual meting. 

The present program also helped us initiate the development of what we believe will be a 
fruitful approach to predicting the important influence of deposition mechanism on resulting 
deposit density and microstructure, for multiphase (eg. granular particulate) deposits. We 
introduced a Monte-Carlo (stochastic simulation) approach to describe particle incorporation events 
at the ’base' of the mass transfer boundary layers (Tassopoulos and Rosner,1991; Tassopoulos, 
O'Brien and Rosner, 1989). In this way, we were able to predict particulate deposit microstructure 
for each condition of particle arrival (mechanism of arrival, particle size distribution, particle 
orientation distribution, particle sticking behavior , etc.). This work, which was made possible by 
our access to the NASA-Lewis Supercomputer, could soon be extended to describe the 
microstructure of deposits formed from vapor mixtures containing suspended particles, say, by 
"PPCVD", setting the stage for analogous treatments of CVD-coating microstructures using 
molecular level simulations. 

2.6 Breakdown of Validity of "No-Slip" Boundary Conditions and Navier-Stokes 

Gas Dynamics in CVD and PVD Applications 

Based, in part, on work carried out under this NASA Grant (and NAG3-898) we have 
also identified the importance of hitherto neglected gas-kinetic phenomena — eg., thermal gas 
"creep" along non-isothermal surfaces and "thermal stress convection" in PVD and CVD 
applications (Rosner, 1989). We demonstrated that these "non-classical" transport phenomena, 
predicted by "higher-order" gas kinetic theories but now known to be important in highly non- 
isothermal low Reynolds number flows, would be especially noticeable (and amenable to study) 
under microgravity conditions. 


3. ADMINISTRATIVE INFORMATION 

The following sub-sections summarize some pertinent 'non-technical' facets of the 
abovementioned Yale HTCRE Lab/NAS A-Lewis research program: 

3.1 Personnel 

The present results (Sections 2 and 5) are due to the contributions of the individuals listed 
in Table 3.1-1, which also indicates the role of each researcher and the relevant time interval of 
their main activity. It will be noted that, in addition to the results themselves, this program has 
simultaneously contributed to the research training (Section 3.2) of a number of students and 
recent PhDs, who will now be in an excellent position to make future contributions to CVD 
technologies oriented toward ceramic coatings. 


- 6 - 



Table 3.1-1 Summary of Research Participants a on NASAGrant : NAG 3-8834 

Research on Chemical Vapor Deposition Processes for Advanced 

Ceramic Coatings 


Name 

Statu s2 

Date(s) 

PrinciDal Research Activity^ 

Castillo, J.L. 

VS 

Summers '88, '89 Asymp theory of vapor phase ignition 

Collins, J. 

GRA 

'88-92 

CVD of ceramic coatings 

Garcia- Ybarra, P. 

VS 

Summers '88/89 Kinetic theory of Soret coeff. 

Kho,T. 

GRA 

92 

chemical vapor infiltration (coatings) 

Rosner, D.E. 

PI 

'88-'92-'93 

program direction-dep.theory/exp 

Tandon, P. 

GRA 

91-92- 

transport phenomena in BLs and CDFs 

Tassopoulos, M. 

GRA 

'88-91 

deposit microstructure simulations 


a PDRA=Post-doctoral Research Asst GRA= Graduate Research Assistant 

PI = Principal Investigator VS = Visiting Scholar 

b See Section 5 for specific references cited in text (Section 2); 

3.2 Presentations and Research Training 

Apart from the publications itemized in Section 5 , our results have also been presented at 
annual or topical conferences of the following professional organizations: 

American Assoc. Aerosol Research (6/90) 

ASME-Engineering Foundation (3/14/91) 

AIChE ( 1 1/9/89, 1 jaw 6/90, 1 1/3/92) 

CVD Xl-Seattle (10/90) 

CVD XII-Hawaii (to occur 5/93) 

Electrochemical Soc., Symposium on Lamp Chemistry (5/88) 

Int. Heat Transfer Conference (Jerusalem) (8/89) 

Materials Research Soc. (1 1/29/89, 12/4/91 
Physicochemical Hydrodynamics (PCH) at MIT 6/27/89 


In addition, during the period: 3/1/88-11/1/92, the PI presented seminars at the following 
Universities: 


Athens (NTU) 10/18/91 
Brown 9/22/92 
Christ Church . (10/91) 

Leeds (5/29/92) 

Limoges (10/10/91) 

Manchester Inst.Sci.Tech. (5/28/92) 
Notre Dame (10/27/92) 

Penn State (7/28/92) 

Sydney (10/11/88) 

Technion (Haifa) 10/27/91) 

Trondheim (9/17/91) 

Waterloo (7/11/91) 

In all, a total of some 32 external talks were 

research program. 


Auckland (1 1/88) 

Brisbane 10/91 
KTH-Stockholm (9/91) 

J. Hopkins U. 

Newcastle (Australia) (11/88) 

Oslo (9/18/91) 

Princeton (9/90) 

Surrey (9/91) 

Toulouse (9/91) 

UNED-Madrid ('89) 

Wisconsin (4/91) 

given based, in part, on the results of this NASA 


This program has involved the PhD dissertation research of three Yale graduate students (J. 
Collins , M. Tassopoulos and P. Tandon; cf. Table 3.1-1). Indeed, J. Collins, who was supported 
for 4 years on a related NASA Training Grant, is expected to complete his Yale U. PhD degree 
requirements in Fall '93. 

Moreover, a Yale Undergraduate ChE major, Jim Tsang (Yale BChE ’93), also completed a 
summer (91) research project on this CVD reactor. 
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3.3 Interaction With US Industry 

The research summarized here was supported by NASA under Grant NAG 3-884 between 
3/1/88 and 11/1/92. During this same period the Yale HTCRE Laboratory has also been the 
beneficiary of continuing smaller grants from U.S. industrial corporations, including Union 
Carbide, SCM-Chemicals, GE-Schenectady, DuPont, and Shell, as well as the feedback and 
advice of principal scientists/engineers from each of these corporations and Combustion 
Engineering-ABB and AVCO-Textron-Lycoming. We appreciate this level of collaboration, and 
expect that it will accelerate inevitable applications of our results in areas relevant to their 
technological objectives in the area of high temperature coatings. 

3.4 Impact/Applications 

It has been particularly gratifying to see direct applications of some of this generic NASA- 
supported CVD research in more applications-oriented investigations reported in recent years. 
Indeed, the writer would appreciate it if further examples known to the reader can be brought to his 
attention. 

In the area of multicomponent vapor deposition in combustion systems additional 
applications of our predictive methods (for "chemically frozen" (Rosner et.al., 1979) and LTCE 
multicomponent laminar boundary layers) continue to be made by British Coal Corporation- Power 
Generation Branch (I. Fantom, contact) in connection with their topping cycles which run gas 
turbines on the products of fluidized bed coal combustors/gasifiers. The writer has also proposed 
applications of our computational methods to Aerojet Propulsion Div. ( D.M. Jassowski ) in 
connection with predicting the chemical stability of iridium rocket nozzle coatings. 

Except for its application to particle-modified CVD systems (eg., Komiyama, 1987 , and 
Hurt and Allendorf, 1991 )particle deposition (PD) was not emphasized in the present NASA 
program. However, for completeness we should mention that explicit PD-examples are provided 
in ongoing work at MIT (Walsh et.al. 1992), and Sandia CRF, both groups having incorporated 
our rational correlation of inertial particle impaction (e.g. a cylinder in cross-flow) in terms of an 
effective Stokes number. This PI was also pleased to confirm recent applications of our deposition 
research (on the correlation of inertial impaction by cylinders in crossflow) by the National 
Engineering Laboratory (NEL) of Glasgow Scotland (Contact: Dr. Andrew Jenkins). NEL is 
apparently collaborating with Marchwood Labs-CEGB on developing mass-transfer prediction 
methods applicable to waste-heat recovery systems in incinerators, as well as pulverized coal-fired 
boilers. These applications are somewhat similar to those reported by the Combustion Lab R&D 
group at MIT and Penn State U. Explicit use of our studies of self-regulated "capture" of incident 
impacting particles (Rosner and Nagaragan, 1987) is also being made in current work on impact 
separators and ceramic heat exchangers for coal-fired turbine systems in high performance 
stationary power plants. Related potential applications arise in connection with "candle filters" used 
to remove fines (sorbent particles,...) upstream of the turbines. A useful summary of work in 
these interrelated areas (Solar Turbines, Textron Defense Systems, Hague International,...) was 
presented at the Engineering Foundation Conference Inorganic Transformations and Ash 
Deposition During Combustion. , the proceedings of which appeared during 1992. 

* * * * 

Clearly, fruitful opportunities for the application of our recent inorganic vapor and/or 
particle deposition research now exist in many of the programs currently supported by NASA, as 
well as US industrial R&D. 
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4. CONCLUSIONS 


Because of (i) the need for tractable predictive methods, instructive prototypical results, and 
rational correlation formulae suited to the needs of the CVD-coatings community, and (ii) our 
previous theoretical and experimental progress in the areas of vapor deposition (VD), particle 
deposition (PD), and their interactions , in March 1988 we embarked on a theoretical progam 
embracing each of the fundamental topics: 

* Role of High Activation Energy Homogeneous Reactions in Limiting 
CVD-Rates and Deposit Quality for Heated Surfaces 

* Prediction of Transport-Shifted CVD-Phase Diagrams; Theory of Deposit 
Composition 

* Role of Vapor Phase Chemistry in Controlling CVD-Rates (LTCE vs CF) 

* Scavenging of Vapors by Suspended Particulate Matter 


* Morphological Stability of Thermophoretically Growing Deposits; 
Transport and Kinetic Phenomena Governing CVD-Coating Surface 
Topography and Deposit Microstructure 

and 

* Breakdown of "No-Slip" Boundary Condition and Navier-Stokes Gas 

Dynamics in CVD and PVD Applications 


Our contributions to the CVD-field in several of these areas (Sections 2, 5.2) were closely 
coupled with corresponding CVD-measurements and calculations at NASA-Lewis Research 
Center. Thus, new concepts and capabilities developed in this program guided the direction of 
closely-related CVD-investigations at NASA-Lewis Labs. This symbiosis benefitted both 
programs— eg., led to respected (and sometimes unforeseen) developments in the area of alkali 
sulfate CVD, some of which are already documented in the References (see, eg. Luthra,1989). In 
each of the specific tasks, the most relevant background literature or precedent is cited (Section 
5.1) . 


The cumulative results of this NASA-supported program, along with complementary 
studies at NASA-Lewis Labs, are expected to constitute a significant body of advances, allowing, 
for the first time, many rational predictions/optimizations for CVD-equipment. The emphasis of the 
fundamental studies described above was on obtaining and communicating an understanding of 
physicochemical phenomena occurring in CVD reactors. In each case we pointed out and explored 
the engineering implications of our Findings in those areas of refractory coating synthesis vital to 
the international competitiveness of the U.S. combustion turbine/chemical rocket industry. 

Based in part on the fruitfulness of the CVD program summarized in the present Final 
Technical Report (and the publications listed in Section 5 below) this laboratory is currently 
exploring opportunities for future collaboration with NASA-Lewis Lab research staff-especially in 
the area of combustion phenomena oriented toward materials synthesis. 
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CVD Chemical vapor deposition CFBL Chemically Frozen BL 
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TTIP Titanium tetra-isopropoxide PD Particle Deposition 

VD Vapor Deposition 
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A NONEQUILIBRIUM THEORY OF SURFACE DEPOSITION 
FROM PARTICLE-LADEN, DILUTE CONDENSIBLE 
VAPOR-CONTAINING LAMINAR BOUNDARY LAYERS 

J. L. CASTiLLot and D. E. RosnerJ 

High Temperature Chemical Reaction Engineering Laboratory, Department of Chemical Engineering, 
Yale University, New Haven, CT 06520, U.S.A. 


( Received 1 January 1986; in revised form 27 June 1987) 

Abstract— The deposition rate of a condensible substance from, say, flowing combustion products to 
“cold" solid surfaces can be strongly influenced by the simultaneous presence of a particulate aerosol since 
the particles can: (a) “scavenge" vapor, thereby influencing the vapor deposition rate; and (b) thermo- 
phoretically drift to the cold surface, carrying their inventory of scavenged condensate. A rational, yet 
quite tractable thermophysical model of these nonequilibrium processes is developed here for high 
Reynolds number laminar stagnation region boundary layer flow, and implemented to the point of 
calculating and displaying the effects of mainstream particle loading, vapor loading and particle size on 
the deposition rate of condensible material at surface temperatures well below the vapor dew point. 
Despite the complexity of this multiphase flow situation, our theoretical model is quite general and is cast 
in terms of dimensionless parameters which dictate the importance of vapor-phase scavenging and particle 
thermophoresis, as well as the Kelvin (surface tension) effect in modifying the submicron particle 
(free-molecule) growth law. Illustrative numerical results are displayed for the deposition of alkali 
sulfate-like vapor from the combustion products of hydrocarbon fuel (or coal) with air, including the 
interesting “structure" of such nonequilibrium multiphase boundary layers. As a useful by-product, our 
results reveal which combinations of particle-phase parameters cause (a) previous “uncoupled" 
vapor/particle deposition rates to be approximately valid as well as (b) recent local vapor/condensate 
equilibrium limit results to be sufficiently accurate. We conclude with an outline of straightforward 
extensions of the present theory to include such factors as: (i) a nonuniform (polydispersed) particle size 
mainstream aerosol; and (ii) size-dependent particle thermophoretic diffusivity; which are likely to be 
important in current or future engineering applications. 


1 INTRODUCTION — MOTIVATION 

The performance of engineering equipment exposed to high-temperature gases containing both 
suspended particulate matter (“dust”) and condensible vapors is determined, in part, by deposition 
phenomena, which can lead to fouling and/or corrosion (e.g, Rosner & Atkins 1983; Rothman 
1985). Even relatively recent work in the areas of submicron particle deposition (e.g. Goren 1977; 
Walker et al. 1979; Gokoglu & Rosner 1985; Rosner & Kim 1984) and vapor deposition (Rosner 
et al. 1979) has, for simplicity, imagined that these processes are separable, i.e. noninteractive. 
However, very recently we have developed a quantitative theoretical explanation of experimentally 
observed deposition rate reductions experienced by highly cooled surfaces due to near -equilibrium 
condensation processes occurring within the gas near the deposition surface (Castillo & Rosner 
1 988a— c; Liang et al. 1988). Actually, in most engineering applications it is not immediately evident 
whether vapor/condensate equilibrium (VCE) can be achieved within such two-phase boundary 
layers (BLs). A more comprehensive theory of such situations would not only yield explicit 
dimensionless criteria for “nearness to VCE”, but would also enable the calculation of fully 
nonequilibrium vapor/particle transport situations — including the important effect of suspended 
aerosol particles on the rate at which condensible material arrives at the cooled surface (Rosner 
& Liang 1988). Our goal here is to develop a rational, but quite general and tractable mathematical 
model for accomplishing these seemingly ambitious goals. 

1.1. Outline of the Present Paper 

In this paper we present a tractable nonequilibrium theory of the deposition on “cold” surfaces 
of unary dilute vapors in high Reynolds number two-dimensional stagnation flow configurations. 

tOn leave of absence from Department of Fundamental Physics, U.N.E.D., Madrid. Spain. 
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BOUNDARY LAYER COAGULATION EFFECTS ON 
THE SIZE DISTRIBUTION OF THERMOPHORETICALLY 
DEPOSITED PARTICLES 1 

H. M. PARK' and DANIEL E. ROSNER' 

Depariment of Chemical Engineering, Yale University, New Haven, CT 06520, U.S.A. 

(Accepted in reused form !0 March 1989) 

Abstract In highly loaded aerosol flow systems of industrial interest, to what extent is the particle size 
distribution (PSD) in a deposit altered by inevitable coagulation events within the thin two-phase “boundary 
layer” adjacent to the collecting surface? As a timely example, under typical optical waveguide preform 
deposition conditions (which include thermophorctically dominated submicron particle deposition at high 
particle mass loadings), we predict that, while Brownian diffusion plays a negligible role in determining the 
total particle mass deposition rate, Brownian coagulation (along with gas shear) plays an important role in 
determining the size distribution of the depositing particles of doped silica. Since the size spectrum of the 
depositing particles, here assumed to be spherical, will strongly influence all thermophysical and structural 
properties of the porous deposited glass, our mathematical model/illustrative calculations can be used to 
circumvent the difficult problem of optically probing the thin laminar thermal boundary layer to infer this 
size spectrum information. Moreover, an understanding of “PSD shifts” induced by coagulation in the 
immediate vicinity of sampling probes/surfaces (including “thermophoretic sampling” recently applied in 
organic soot research) can be used to correct electron-microscope-inferred particle size/shape distributions 
for the systematic effects of coagulation within such boundary layers in systems with higher particle mass 
loadings. 


INTRODUCTION 

In most previous particle deposition rate theories [see, 
for example, Friedlander (1977) and Rosner (1989) 
(who explicitly considered convection/Browmian dif- 
fusion); and Goren ( 1 977), Rosner and Park ( 1 988) and 
Park and Rosner (1986, 1987) (who considered 
thermophoretically driven submicron aerosol particle 
transport through nonisothermal combustion prod- 
uct gases)] it has been assumed that the suspended 
particles were spherical and remained constant in 
diameter up to the moment of deposition. However, as 
is well known, aerosols are unstable with respect to 
coagulation, and, especially in highly loaded aerosol 
systems, coagulation may exert an important influence 
on system or product performance. In previous in- 
vestigations of thermophoretic deposition rates under 
conditions of high particle mass loading we could 
neglect particle coagulation phenomena by exploiting 
the fact that, since the thermophoretic diffusivity, 
x T D p remains nearly constant with respect to the 
particle Knudsen number for submicron particles in 
atmospheric pressure combustion products (Rosner, 
1 980), particle coagulation does not affect the resulting 


f Originally submitted May 1987 to Chemical Engineering 
Science. Based, in part, on H. M. Park’s PhD dissertation, 
Department of Chemical Engineering, Yale University 
(1987). 
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aerosol deposition rate (Rosner and Park, 1988). How - 
ever, in many industrial applications, including the 
manufacture of “preforms" for optical waveguides, the 
size distribution of the depositing particles will exert a 
strong influence on the porosity and other important 
microstructural properties of the deposit 
(Tassopoulos et a/., 1989), thereby affecting all sub- 
sequent processing (e.g. drying, spinning etc.) (Keiser. 
1983). As a useful first step in the prediction of size 
distribution “shifts" for thermophoretically deposited 
particles we present here a simple mathematical model 
that describes particle coagulation behavior across 
nonisothermal highly-loaded “dusty-gas" laminar 
boundary layers (BLs). 

THEORETICAL FORMULATION 
Overall BL flow structure 

If we neglect particle inertial drift [see Rosner and 
Park (1988, Appendix)], the Eulerian BL flow equations 
for highly particle mass loaded systems treated as a 
continuum are [see, for example, Fernandez de la 
Mora and Rosner (1982) and Rosner and Park (1988)] 
a mixture mass balance, x-momentum balance, mix- 
ture energy balance and total particle mass (irrespec- 
tive of particle size) balance. In their self-similar forms 
these equations, combined with an appropriate mix- 
ture equation of state {p = [pMJ(RT)']{ \ -cu p )“ l } and 
boundary conditions, were derived and solved nu- 
merically by Rosner and Park (1988) to obtain the 
overall boundary layer structure and wall transfer 
rates of x-momentum, mass and energy, under para- 
metric conditions relevant to optical waveguide pro- 
cessing. It was pointed out [Rosner and Park (1988, 
Appendix 5)] that, while the explicitly considered 
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Side-wall gas “creep” and “thermal stress convection” In microgravity 
experiments on film growth by vapor transport 

Daniel E. Rosner 

Department of Chemical Engineering, High Temperature Chemical Reaction Engineering Laboratory \ 

Yale University, New Haven, Connecticut 06520 

(Received 13 July 1989; accepted 15 August 1989) 

While “no-slip” boundary conditions and the Navier-Stokes equations of continuum fluid 
mechanics have served the (moderate-to-“high” pressure) vapor transport community well 
until now, it is pointed out that transport conditions within highly non isothermal ampoules 
(e.g., those used to grow organic solid thin films in microgravity experiments) are such that 
the nonisothermal side walls “drive” (rather than damp) the dominant convective flow, and 
the familiar Stokes-Fourier-Fick laws governing the molecular fluxes of momentum, energy, 
and (species) mass in the “continuum” field equations will often prove to be inadequate, even 
at Knudsen numbers as small as 10“ 3 . The implications of these interesting gas kinetic 
phenomena under microgravity conditions, and even under Earth-bound experimental 
conditions, are outlined here, along with a tractable approach to their systematic treatment. 


It is well known that even in outwardly simple, closed 
ampoules, rather complex vapor motion can occur as a result 
of the effects of a body force (usually gravity 1 ) and/or in- 
duced by the net mass transfer itself (the so-called Stefan 
flow). 2 ' Using familiar “no-slip” boundary conditions, 
these have been analyzed in simple situations (ampoule ge- 
ometry and orientation with respect to gravity, small frac- 
tional change in mixture density, equal vapor diffusivities for 
the transport of momentum, energy and mass, no side-wall 
effects on energy and mass transfer, etc.), with a typical set 
of results reported by Miller. 4 We recently completed a pre- 
liminary theoretical study whose goal was to converge on the 
simplest thermophysical model capable of accurately de- 
scribing film growth conditions in sealed nonisothermal am- 
poules used to grow organic thin films under microgravity 
conditions (see, e.g., Refs. 5 and 6). In the course of this 
study, we discovered* that, whereas the above-mentioned 
Stefan- and buoyancy-driven flows are likely to be negligible 
under the typical microgravity conditions employed, instead 
several previously neglected “gas kinetic” phenomena (con- 
sequences of the Boltzmann equation of low-density gasdy- 
namics) 7 become rather important “sources” of convection. 
These less familiar phenomena are discussed sequentially be- 
low. The importance of non-Fickian vapor diffusion via the 
Soret H (thermal mass diffusion) effect is also briefly consid- 
ered. We conclude with the future implications of these in- 
teresting departures from the conventional Navier-Stokes- 
Fourier-Fick convection-diffusion laws, with emphasis on 
developing a tractable approach to the systematic inclusion 


of these gas kinetic effects in the general field of physical 
vapor transport (PVT) and chemical vapor transport 
(CVT). 

To the best of our knowledge, in all previous analyses of 
vapor transport (“chemical” or “physical”) in nonisother- 
mal ampoules, it has been explicitly or implicitly assumed 
that the ( motionless ) ampoule side walls will bring any local 
gas tangential motion to rest — a familiar hypothesis known 
as the “no-slip” condition. From this “classical” perspective 
the side walls “damp” any gas motions established as a result 
of gravity and/or the net vapor transfer process itself. How- 
ever, the quantitative estimates below reveal that, in fact, the 
nonisothermal side walls will act like a “pump,” producing a 
side-wall gas “creep,” which would be present (indeed 
dominant ) in the complete absence of gravity and/ or Stefan 
flow! For this reason, we outline here the essential features of 
such nonisothermal side-wall gas “creep” and call attention 
to its implications for future vapor transport microgravity 
experiments and ancillary modeling efforts. Even in ground- 
based experiments on ampoules at the same pressure level, 
we estimate that side-wall-induced gas “creep” velocities 
can be within one order of magnitude of those associated 
with conventional buoyancy-induced gas velocities — con- 
vection known to strongly influence “solute” vapor trans- 
port rates, and film quality (uniformity, orientation, etc ). 

Just over one century ago, and in the last year of his life, 
Maxwell pointed out 9 that (unless gas/solid encounters 
were completely “specular,” i.e., with no tangential momen- 
tum exchange) a nonuniform temperature surface would in- 
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PREDICTION AND USE OF TRANSPORT-SHIFTED CVD PHASE 
DIAGRAMS: TITANIUM DIBORIDE CVD FROM TICI4/BCI3 MIXTURES* 


D.E. Rosner* and J. Collins b 

High Temperature Chemical Reaction Engineering Laboratory 
Yale University Dept, of Chemical Engineering. 

New Haven CT 06520-2159. USA 


ABSTRACT 

An efficient method for predicting mass "transport-shifted" CVD 
phase diagrams Is described. Prototypical results are presented 
revealing Important systematic shifts in deposit phase/stoichiometry 
due to chemical element segregation in the mass transfer controlled 
CVD of TiB 2 . Our convective-diffusive mass transfer model accounts 
for both multicomponent Fick and Soret transport of dilute species 
(18 for TIB 2 deposition) across the gaseous boundary layer, with local 
thermochemical equilibrium imposed at the deposition surface. 
Selected results of this analysis, which also yields loci of constant rates 
of deposition or etching, are compared with available experimental 
data within the predicted domain of self-consistency (i.e. under 
conditions such that neither homogeneous reactions nor 
heterogeneous kinetic barriers are expected to play a significant role). 


1 . INTRODUCTION 

In chemical vapor deposition (CVD) systems the elemental 
composition established at the gas/solid interface is generally different 
from the feed stream elemental composition. Therefore, conventional 
CVD phase diagrams based on equilibrium calculations using feed 
stream element fractions cannot be expected to yield the equilibrium 
state corresponding to the temperature, pressure, and element 
fractions actually established at the gas/solid interface, where the 
decisive CVD chemical reactions occur. Indeed, mass transport 
induced shifts in elemental composition between the feed stream and 
the gas/solid interface in CVD systems have often been suggested as a 
qualitative explanation for experimentally observed shifts in CVD 
deposit phase boundaries, as in the S1B X deposition studies by Spear 
and Dirkx (1.2). However, we are unaware of any systematic treatment 
of the required multicomponent mass transfer problem. 


* Professor of Chemical Engineering and Director. Yale HTCRE Laboratory 

b Graduate Research Assistant and NASA GSRP Fellow, Yale HTCRE Laboratory 
Research Supported by NASA Lewis Research Center Grants NTG-50377, NAG 3-884. 

• Submitted 5/4/90 for presentation at the Eleventh International Conference on 
Chemical Vapor Deposition ((CVD XI). Seattle WA, October 14-19. 1990. 
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Morphological instability of a thermophoretically growing deposit 

Jose L. Castillo, Pedro L. Garcia-Ybarra 
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and 

Daniel E. Rosner 
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C'onneein ut 0ft520-2l59YS, USA 

Received 15 May 1990; manuscript received in final form 15 April 1991 


The stability of the planar interface of a structureless solid growing from a depositing component dilute in a carrier fluid is 
studied when the main solute transport mechanism is thermal (Soret) diffusion. A linear stability analysis, earned out in the limit of 
low growth Peclet number, leads to a dispersion relation which shows that the planar front is unstable either when the thermal 
diffusion factor of the condensing component is positive and the latent heat release is small or when the thermal diffusion factor is 
negative and the solid grows over a thermally-insulating substrate. Furthermore, the influence c|f interfacial energy effects and 
constitutional supersaturation in the vicinity of the moving interface is analyzed in the limit of very^£2& Schmidt numbers (small 
solute Fickian diffusion). The analysis is relevant to physical vapor deposition of very massive species on cold surfaces, as in recent 
experiments of organic solid film growth under microgravity conditions. 


1. Introduction 

The growth of solid deposits from the vapor 
phase is a technique widely used for the fabrica- 
tion of thin layers of metals, insulators, and semi- 
conductor materials; see, for instance, the review 
papers by Stringfellow [1] and Bryant [2], Even 
though crystal growth from melts may be a more 
“efficient" process (usually giving growth rates 
several orders of magnitude larger than vapor 
phase growth), the latter is often preferred be- 
cause of film quality considerations. In particular, 
new pharmacological and electro-optical applica- 
tions demand the production of crystalline films 
of organic and organometallic macromolecules. 
The optimal growth of these materials from the 
vapor phase is now- under intense experimental 


investigation internationally, even in microgravity 
environments to avoid crystal imperfections in- 
duced by buoyancy-driven instabilities, see Wal- 
ter [3]. 

Especially since the pioneering work of Mullins 
and Sekerka [4,5] and the review by Langer [6], it 
is well known that planar crystal growth may be 
morphologically unstable during the solidification 
of a supercooled liquid (i.e., in the absence of 
capillarity, a forward-facing bulge in the interface 
steepens the thermal gradient in the fluid ahead 
of it, so that latent heat is more rapidly removed 
from the surface and the bulge amplifies unsta- 
bly), and in the similar case of isothermal solidifi- 
cation of a liquid mixture (respectively, direc- 
tional solidification of a binary liquid) due to the 
constitutional supersaturation (respectively, const i- 
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ONSET CONDITIONS FOR GAS PHASE REACTIONS AND PARTICLE 
NUCLEATION/GROWTH IN CVD BOUNDARY LAYERS. 

D.E. Rosner® , J. Collins' 5 and J.L. Castillo 0 
High Temperature Chemical Reaction Engineering Laboratory 
Yale University, Department of Chemical Engineering 
New Haven CT 06520-2159, USA 

Literature CVD-rate data, and our recent experiments on Ti02(s) 
film growth from titanium (IV) tetra-isopropoxide (TTIP) vapor 
using a well-defined impinging jet reactor, reveal that the onset 
of vapor phase reactions near a hot deposition surface can lead to 
sharp reductions in CVD-rate and alterations in deposit 
microstructure. These observations have motivated our 
development of a thin chemical sublayer (CSL-) theory for 
predicting (using simple formulae) the interplay of 
heterogeneous kinetics, homogeneous kinetics and (Fick-, Soret-, 
convective-) transport phenomena in CVD reactors. CSL theory" 
only briefly outlined here, can be used to interrelate different 
CVD reactors /conditions and thereby guide the selection of 
reactor configurations/conditions that will lead to the 
maximum CVD-rate prior to the onset of vapor phase reactions 
and/or external transport limitations. Comparisons with our 
present Ti02(s) ( via TTIP) data are used to sutyjest fruitful 
extensions of this work. 

1. INTRODUCTION 

To interpret our recent experimental observations on the CVD rate of 
T 1 O 2 films in an impinging jet reactor (1,5) and guide the design /operation of 
future CVD reactors, a rational asymptotic theory has been developed for the 
onset of homogeneous chemical reactions in the thermal boundary layer (BL) 
adjacent to a hot deposition surface, with emphasis on their effect on 
deposition rates and film quality (1,2). The analysis is tailored to systems in 
which a) gas phase reactions are confined to a thin chemical sublayer (CSL) 
embedded inside the thermal BL and b) the homogeneous reaction(s) of 
interest can be represented by a single high activation energy (Ehom/(RT e )» 1) 
Arrhenius-type rate expression. In this limit one finds (2): 



* Professor of Chemical Engineering and Director, HTCRE Laboratoiy HTCRE Ms #195 
b Graduate Research Assistant, Yale HTCRE Laboratory * 
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